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Abstract 



Competition among particle evaporation, temperature gradient and flow is 
"q ■ investigated in a phenomenological manner, based on a simultaneous anal- 

ysis of quantum statistical correlations and momentum distributions for a 
non-relativistic, spherically symmetric, three-dimensionally expanding, finite 
source. The parameters of the model emission function are constrained by fits 
to neutron and proton momentum distributions and correlation functions in 
■ intermediate energy heavy-ion collisions. The temperature gradient is related 

to the momentum dependence of the radius parameters of the two-particle 
correlation function, as well as to the momentum-dependent temperature pa- 
rameter of the single particle spectrum, while a long duration of particle evap- 
oration is found to be responsible for the low relative momentum behavior of 
the two-particle correlations. 
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I. INTRODUCTION 



It has been shown recently that for non-relativistic, three-dimensionally expanding sys- 
tems the quantum statistical correlations measure only part of the particle source and that 
the effective temperature of the momentum distribution is obtained as a combination of the 
freeze-out temperature and of the "geometrical temperature" , a term due to expansion and 
to the finite geometrical sizes 0. The expansion makes the effective radius parameters of 
the two-particle correlation functions smaller than the geometrical size of the source even 
for non-relativistic expanding systems 

Here we extend our study of non-relativistic expanding particle sources to the case of 
a temperature gradient inside the source. Such an extension of the model of Ref. QXJ is 
motivated by data on momentum distribution of protons in heavy ion reactions in the non- 
relativistic energy domain, which show deviations from the purely exponential spectrum 
|J|3|], although the neutron spectrum is well approximated by a Boltzmann distribution. 

The purpose of our paper is to investigate the features of non-relativistic heavy ion 
collisions at a phenomenological level, parameterizing the emission function assuming local 
thermal equilibrium, without attempting to describe the microscopic mechanisms which lead 
to such a particle emission pattern. However, we include phenomenologically general features 
like evaporation, cooling, temperature gradient and flow as well as strong and Coulomb final 
state interactions for protons and neutrons, and test the resulting model on data at E = 30 
MeV/nucleon. Analytic approximations are also derived in order to gain improved insight 
into the influences of the competing effects. These analytic results are not directly applicable 
to nucleon correlations. However, they may be applicable to study pion correlation functions 
and spectra in such reactions where a spherical non-relativistic source may be assumed. 

At high energies a relativistic analogue model j|] has been successful in describing the 
available data ||. However there are important qualitative differences between relativistic 
heavy ion collisions at CERN SPS and those at non-relativistic energies from the point of 
view of particle sources. Low and intermediate energy reactions may create a very long- 
lived, evaporative source, with characteristic lifetimes of a few 100 fm/c, in contrast to the 
relatively short-lived systems of lifetimes not larger than a few 10 fm/c at CERN SPS. During 
such long evaporation times, cooling of the source is unavoidable and has to be included 
into the model, while the freeze-out at a constant temperature is a much better description 
at higher energies. Furthermore, in the non-relativistic heavy ion collisions mostly protons 
and neutrons are emitted and they have much stronger final state interactions than the 
pions dominating the final state at ultra-relativistic energies (for recent reviews on nucleon 
interferometry, see for example Refs. @-|8[]). 

We would like to emphasize that it is not our purpose to re- invent a detailed microscopical 
description of non-relativistic nuclear collisions. We are trying to develop a framework to 
describe the particle emission patterns, i.e., a simultaneous description of the invariant 
momentum distribution and the quantum statistical correlation function, from intermediate 
energies up to the highest energies achievable. In this paper we concentrate on the low 
part of the intermediate energy range. For example, sophisticated microscopical transport 
descriptions f9j, such as the BUU (Boltzmann-Uehling-Uhlenbeck) and the QMD (Quantum 
Molecular Dynamics) models, are well-known and believed to provide a reasonable picture of 
proton emission in central heavy ion collisions from a few tenths up to hundreds of MeV per 
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nucleoli. However, in Ref. [10|, the BUU model predicts too large correlations and under- 
predicts the number of protons emitted with low energies, for the reaction 36 Ar + 45 Sc at 
E = 120 and 160 MeV/nucleon. This indicates that the simultaneous description of two- 
particle correlations and single-particle spectra is a rather difficult task. For energies below 
a few tens of MeV per nucleon, where long-lived evaporative particle emission is expected to 
dominate, the measured two-proton correlation functions were found to be consistent with 
compound- nucleus model predictions |]TTJ; however, a simultaneous analysis of proton and 
neutron single particle spectra and two-particle correlation has not yet been performed to 
the best of our knowledge. 



The basic model, including a temperature gradient, is presented in section |11 A|. Section 



I [B] contains an extension of the model to long emission times. Analytic approximations for 
momentum distributions and correlation functions are derived in section [Til], and applied in 
section [TV] to situations where final state interactions can be ignored. Section |V] contains 
numerical applications to neutron and proton interferometry, where final state interactions 
are important. Finally, our results are summarized in section [VI], while appendix [A] contains 
additional material on analytical approximations for the momentum distribution. 



II. THE MODEL 

The model presented in this section is based on the work of Ref. jj]]. For clarity we briefly 
summarize the formalism along the lines of Ref. |1| , together with a thorough description of 
the extensions. 



A. Emission Function with a Temperature Gradient 

The emission is characterized by the emission function or source function S(x,p) which 
is the probability that a particle is produced at a given x — (t, r) — (t,r x ,r y ,r z ) point 
in space-time with the four-momentum p = (E,p) = (E,p x ,p y ,p z ), where the particle is 
on mass shell, m 2 = E 2 — p 2 . The quantum- mechanical analogy to the classical emission 
function is the time-derivative of the non-relativistic Wigner-function [0| which is analogous 



to the covariant Wigner-function of Ref. [|l3j. In this work the time-derivative of the Wigner- 
function will be approximated by classical emission functions. 

In terms of the emission function both the invariant momentum distribution (IMD) 
and the quantum statistical correlation function (QSCF) are prescribed. The Fourier- 
transformed emission function is introduced as an auxiliary quantity 

S(Ak; K) = J d 4 x S(x; K) exp( iAk • x) , (1) 

where 

Ak = p l -p 2l K = — - — (2) 

and Ak-x stands for the four-product of the four- vectors. Then the momentum distribution 
of number of the emitted particles, iV\ (p) is given by 
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^Px Uj i J y ^yz 

which is normalized to the mean multiplicity as 

S V N x (p) = (n). (4) 



Note that this normalization condition is different from the one used in Ref. [1] where 
the momentum distribution has been normalized to unity. We also assume that the non- 
relativistic measure of the momentum space should be applied, due to the non-relativistic 
nature of the considered problem. 

In the plane- wave approximation (i.e. neglecting final state interactions), the Bose- 
Einstein or Fermi-Dirac correlation functions are prescribed in terms of our auxiliary function 
as 

where the + sign stands for bosons and the — sign for fermions. This approximation involves 
an off-shell continuation of the on-shell emission functions, the significance of which was 



discussed first in Refs. 0,113 



For central heavy ion collisions at intermediate energies the target and the projectile form 
a collective state which can be described as a non-relativistically expanding fluid within the 
framework of hydrodynamical models. Due to the expansion, the fluid cools and we assume 
that it then undergoes a certain disintegration process. In case the information about 
the initial directions is lost, the final freeze-out stage becomes approximately spherically 
symmetric. 

We assume along the lines of Ref. |I[ that the emission function is characterized by a 
distribution of production points I(r ) and by a distribution of the particle emission times, 
H(t). The correlations between space-time and momentum-space shall be introduced by 
a non-relativistic momentum distribution. We assume that the expanding system is di- 
lute enough when the particles are emitted so that the quantum statistical single-particle 
distribution can be well approximated by a Boltzmann distribution, 

f&P) = exp [ 2mT{x) If (or), (6) 




J(x)=exp • (7) 



Here g is the degeneracy factor, u(x) is the (non-relativistic) flow velocity, the freeze-out 
temperature is denoted by T and n(x) is the chemical potential. 
Thus the emission function is characterized as 

S(x;K) = f(x;K)H(t). (8) 

In order to simplify the results we shall keep only the mean and the width of the source 
distributions, i.e., we shall apply the Gaussian approximations for the distribution functions 
of t and r as follows 
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/(x)=exp( v -jexp^^j, (9) 

H(t) = 7 i_ exp ) . (10) 

v ; (2vrAt 2 ) 1 /2 p ^ 2 At 2 / v ; 



(The expression in Eq. fllOD is not suitable for very long emission times and evaporative 
processes. Such scenarios will be treated and discussed in the next section.) 
In other words, we have the following ansatz for the chemical potential n(x): 



H{x) /i r 
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T(x) T 2R { 
which is analogous to the ones used in Refs. [fLl|l5 



2 ■ 

G 



(11) 



We assume the following form of the local temperature distribution: 

r « = TT^M ■ (12) 

The parameter a controls the gradient of the local temperature at the last interaction point. 
We will in this work only treat the case when the temperature profile is decreasing as a 



function of r, though the expression ( |12|) contains also increasing temperature profiles by 
taking a pure imaginary. Note, that the particular choice of the temperature profile is not 
influencing the leading order, approximate results: namely, we shall apply a saddle-point 
approximation and an expansion of the inverse temperature profile around the r = 0, t = t 
point. Thus any other temperature profile, which leads to similar expansion coefficients, 
shall lead to similar approximate results. 

We assume that the freeze-out temperature at each emission point is much smaller than 
the mass of the particles in this non-relativistic case, 

T(r,t)<m. (13) 

We select a velocity of the 3D expanding matter at space-time point x so that it be 
spherically symmetric and describe an expansion in all three directions with a constant 
gradient. Thus the velocity around the mean freeze-out time to is assumed to have the form 

u(x) = b^, (14) 
which describes a scaling solution of the non-relativistic hydrodynamical equations at the 



mean freeze-out time to fo r | r |<C to an d 6 = 1, Ref. |L6|. Here we have introduced the 
parameter b which controls the amount of flow. For b = we recover the case without flow. 
The results given in Ref. fl[] are re-obtained for the case a = and 6=1. 

In Ref. Jl[] it was implicitly assumed that the duration of the particle emission is short, 
At <^ t , since the flow field, the geometrical radius and the freeze-out temperature were 
assumed not to change significantly during the time interval At, centered on t . The duration 
of the particle emission, At, has thus to satisfy 

A t « mi n(^) (15) 
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in order to warrant the model assumptions. If we have o,5 < 1 then At pa to becomes 
possible since we approach the static fire-ball case. 

When the above approximations for the emission function are valid, the auxiliary function 
can be rewritten as 

S(Ak, K) = H(AE) J d 3 r exp {-ik-r) f(t Q , r; K), (16) 

where H(AE) stands for the Fourier-transformed freeze-out time distribution. Within this 
approximation the freeze-out time distribution determines the energy-difference dependent 
part of the correlation function. 



B. Extension to Long Emission Times 

For low and intermediate energy nuclear reactions the duration time for the particle 
emission can be quite long. Evaporative models (with long duration times) have been fairly 



successful in describing certain aspects of such reactions |17|J18| . The model presented above 
can be extended to also include long emission times by introducing cooling of the source and 
by modifying the time distribution. 

The effect of possible cooling of the source can be caused, e.g., by the expansion and by 
the evaporation of the particles from the source. Although cooling can be consistently in- 



corporated in the model of section [II A| (based on hydrodynamical results) it is very difficult 



to estimate analytically the effects for a long-lived, slowly cooling source, when particles are 
emitted during a large time interval. We therefore instead choose to utilize a phenomeno- 
logical cooling profile, i.e. we assume that the decrease of the local temperature is given 
by 

T(r,t)=T(r)(^y , (17) 

where T(r) is given by Eq. QT^ ) and we use values r = 100 fm/c and a = 1/3, motivated by 
a profile for isentropic three-dimensional expansion of and ideal gas. Note, that this form 
may not be a good approximation for the initial stage of the reaction since such a form 
cannot describe the initial rise of the temperature. However, we numerically found, that a 
long particle evaporation time is needed to get a simultaneous description of particle spectra 
and correlations in intermediate energy heavy ion collisions. Thus, most of the particles are 
emitted much later than the initial stage of the reaction and the approximate profile given 
by the above equation may be suitable to describe the fall of the temperature during most 
of the particle emission. A more detailed temperature profile could be obtained by a fit to 
the local temperature distribution of a microscopical simulation including collective effects 
(re-scattering) and particle evaporation. 

The Gaussian time distribution in Eq. fllQl) does not give a good description for small 
times, t pa 0, when At is large, and thus has to be modified. Closely related to the Gaussian 
approximation is a distorted Gaussian distribution, which vanishes for very small times. We 
assume the form 

H r {t) = 2td exp(- 1 2 d)9(t) , (18) 
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which has the mean emission time 



W = ~ 2 f d (19) 

with the variance 

a 2 (t) = (t 2 )-(t) 2 = l(l-^) . (20) 

This form is motivated by the fact that there must be an initial rise in the particle emission 
followed by a long exponential tail. Such a distribution is for example the gamma distribution 
in t 2 , which contains two parameters. However, for very broad time distributions, the mean 
and the variance are not very sensitive to independent tuning of these two parameters. Thus 
considering the current precision of the data (that is discussed in section |V|), we have fixed the 
small t behavior of the H-p(t) curve to a linear rise, decreasing the number of free parameters 
by one and simplifying the numerical integrations at the same time. Note that the analytic 
expressions above for (t) and <J 2 {t) are valid only when cooling is excluded. When cooling is 
included the mean emission time and the variance will be somewhat modified (see also the 
results and discussion in section |V|). 

The source function S(x,p) with the time distribution, H r (t) in Eq. flIS|), describes an 
expanding source with long duration of the particle emission, including cooling of the source, 
Eq. ([H]). It thus describes the gross features of particle evaporation, which is believed to 
take place at low collision energies. However, as an evaporative model, some features are 
only approximately treated. For example, quantities like flow velocity (u) and geometrical 
source size (Rg) are taken time independent, while in a more rigorous treatment also such 
quantities would vary with the time. Examples of more refined evaporative models can be 
found in Refs. fT7j , [19|j . Here, however, we want to utilize a rather simple model, containing 
few parameters, applicable in a wide energy range (at high energies particle evaporation is 
negligible). We want a model that can describe nucleon emission as well as pion emission, 
and we want the model to be simple enough so that certain analytic results can be derived. 
Thus when the model is applied near its low energy limit (as in section |V|), the extracted 
quantities like flow velocity and source radius, will reflect time averages. 



III. ANALYTIC APPROXIMATIONS 

The integral for the auxiliary function S(Ak, K) in equation ([16]), can be evaluated with 



the help of the saddle-point method, which is described in details in Refs. flq , |20|j21| . At the 
saddle-point, the partial derivatives of the emission function with regard to either r x ,r y or 
r z vanish simultaneously. 

The saddle point, r s is found by assuming that | r s \ /t -C 1. In this case one can expand 
the saddle point equations around r = and one can solve the saddle point equations in a 
linearized problem. Note that for a = the method yields exact results because the non- 
linearity of the saddle-point equations is related to the non-vanishing values of a. Since we 
utilize an expansion for | r | /to <C 1, all shapes of the temperature profile which lead to the 
same second order expansion for small distances as Eq. (O) lead to the same results. 



With the help of the above approximations, we may rewrite the emission function as 
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with 

" = (2tt)3 — VT 
The mean emission point coincides with the saddle-point r s (p), being 



^-T^expf^'l . (22) 



r '< P > = " *> *E>V) + » m + HT t IBV (23) 

where the kinetic energy is denoted by E k (p) = p 2 /(2m). The following result is obtained 
for R*: 

' + Wt^ ( 24 ) 



^ = ^JW ^ 

Within the above approximations, the auxiliary function contains a momentum-dependent 
factor which shall enter the momentum distribution only and another factor which is a 
Fourier-transform of a Gaussian and so easily integrable. The approximations are self- 
consistent if the condition | r s \ /t <C 1 is satisfied. We return to this point later. 

The analytic results for the momentum distribution and the quantum statistical corre- 
lation function are given as 

N 1 (p)=c a (2«R.(p)) exp(^ 2mT{rs{p)) , (26) 

C(K, Ak) = l± exp(-Rl(K)Ak 2 - At 2 AE 2 ) . (27) 

The effects of final state Coulomb and Yukawa interactions on the two-particle relative 
wave-functions are neglected (in section [V] final state interactions are taken into account). 

These expressions are generalizations of the momentum distribution and correlation func- 
tion obtained in Ref. M and reduce to those for a = and b = 1. The static fire-ball corre- 
sponds to the case a = b = 0. A new feature for a ^ is that the radius parameter of the 
QSCF becomes a decreasing function of the momentum and also the effective temperature 
becomes momentum dependent. 

Similarly to the a = and 6=1 case Q, the effective temperature T* shall be determined 
by the maximum of the local temperature, To, and the geometrical temperature defined as 



R, 



The relationship is given by 



T G = T -± . (28) 

Itrn 



1 f 1 - f 

% T + T G T ' K J 

f = (30) 
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FIG. 1. Comparison of the numerically integrated energy distribution of protons (solid line) 
with a Boltzmann distribution of = 23.1 MeV (dashed line) and with analytical results for the 
linearized saddle-point calculation, Eq. ( p6|) as indicated by the dash-dotted line, for the parameter 
values of m = 938.3 MeV/c 2 , R G = 7.6 fm, T = 10.0 MeV, to = 45.0 fm/c, a = 1.0, and b = 0.7. 
The energy distributions, f(Ek) oc Ni(p)/p 2 , are re-scaled so that f(Ek,o) = 1 for E^p = 10 MeV. 
Note the characteristic bending down of both the solid and the dash-dotted line as compared to 
the Boltzmann distribution. 



which is analogous to the case obtained for the slope of the momentum distribution at 



high energies [§,f24]. In the case of no temperature gradient, T* will grow linearly with 
the particle mass. In case of no flow, is independent of mass. For a given mass, the 
effective temperature T* shall be constant only in a limited p interval, which is given by 
p 2 ^ m 2 b 2 /a 2 . Since the present investigation is limited to the non-relativistic p m 
region, it follows that the effective temperature is not noticeably momentum-dependent for 
small temperature gradient satisfying a 2 < b 2 . As soon as the temperature gradient increases 
above the a 2 < b 2 region, a characteristic high momentum suppression may appear in the 
tail of the distribution, see Fig. [I] for example. 

When fitting the model to data, we concentrate on the a 2 > case, since this is the 
case which may result in a suppression at high momentum - a phenomenon observed in the 
momentum distribution of protons in intermediate energy heavy ion reactions [0,0] • 

A short discussion on the validity of the saddle-point approximation is presented in 
appendix together with suggested improvements of the approximation and alternative 
approximation schemes. 
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IV. CORRELATIONS WITHOUT FINAL STATE INTERACTIONS 



In this section the analytic results will be applied to momentum distributions and cor- 
relation functions where final state interactions can be neglected. Applications could for 
example be pion distributions and correlations at BEVALAC energies, though the discus- 
sion in this section is quite general and could be applied also to other situations where final 
state interactions can be ignored. The qualitative discussion in this section also serves as a 
preparation to section [V], where final state interaction is introduced as well. The presentation 
will follow the lines of Ref. JIJ, since the results are very similar. 

The relative momentum, Afc, which appears in the correlation function, is invariant 
under Galilei-transformations, but the energy difference is not invariant even under the non- 
relativistic Galilei transformations. This can be re-formulated so that the specific directional 
dependence becomes more transparent. The energy difference is 

AE = Pl_Pl = (El±Pl).^H = ^ = VK . Ak , (31 ) 

2m \ 2 J m m 

where we have introduced the mean velocity of the pair, vk = K/m. 

Let us define the out direction to be parallel to the mean velocity of the pair, vk, and 
the perp index for the remaining two principal directions, both being perpendicular to the 
out direction. This naming convention corresponds to the one used in high energy heavy ion 



collisions [E21E3 



By this definition, the principal directions, which will be utilized to evaluate the cor- 
relation function, are assigned to the given mean momentum of a given particle pair and 
not that of the fire-ball: by changing the mean momentum we change the out and the perp 
directions as well. In any given frame, the relative momentum can be decomposed as 

Afc = Q ou t + Qperp (32) 

and the correlation function can be rewritten as 

C(Q out , Q perp ) = 1 ± exp (-Rl erp (K)Q 2 perp - R 2 out (K)Q 2 out ) , (33) 

RU(K) = Rl { K) = (34) 

R 2 out (K) = Rl(K) + v 2 K At 2 . (35) 



Thus the perpendicular components shall measure the effective radius R*(K). What is the 
interpretation of this quantity? 

There are two length-scales in the problem: the geometrical length-scale Rq and the 
thermal one Rt, the latter being generated by the flow gradient b/to, the central temperature 
T and the temperature gradient a/t . From the previous equation one can see that the 
effective R*(K), measured by the perp component of the correlation function, is dominated 
by the smaller of the two. In other words, for large and relatively cold three-dimensionally 
expanding systems, not the whole source can be seen by quantum statistical correlations, but 
only a part of the whole system, which is characterized by a thermal length-scale Rt{K). 
Note also that, unlike in the case discussed in Ref. [|1], the radius parameter R*(K) shall 
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FIG. 2. Kinetic energy dependence of the effective radius parameter R*(Ek), Eq. (|2^ - p5|) . 
Solid line indicates R*(Ek) for the parameter values of Fig. dash-dotted line stands for a = 2.0, 
b = 0.1, other parameters being the same as for the solid line. 

not be a constant but it will depend on the momentum of the particles. See Fig. ^| for 
illustration. 

From Eq. (|35| ) it follows that the out component in general shall also be sensitive to 
the duration of the freeze-out time distribution, since it contains a term v^-At 2 . This term 
vanishes in the center of mass frame (c.m.f.) of the particle pair, since AE = \ vk |= in 
c.m.f. In this specific system, 

Q) = Qler P + Qlut - A^ 2 = Qler P + Qlut in c.m.f. of the pair (36) 

is the invariant momentum difference. The correlation function in the considered case be- 
comes symmetric, when evaluated in the cm. of the pair: 

C(Qi) = 1 ± exp(-Rl(K)Qj) in c.m.f. of the pair. (37) 

It is interesting to investigate another limiting case, R~g <C Rt{K). In this case we 
re-obtain the standard results: 

R lerp = R G > R lut = R G + V 2 K At 2 , T* = T , (38) 

i.e., if the thermal length scale is larger than the geometrical size in all three directions, 
the correlation measurement determines the geometrical sizes properly, and the momentum 
distribution will be determined by the freeze-out temperature being just a thermal distribu- 
tion for a static source. In this limiting case, the effective radius parameter as well as the 
effective temperature become independent of the momentum of the emitted particles. 

The non-relativistic model presented in this work has six input parameters: Rq, to, At, 
T, a and b, which determine two measurable functions. This, in turn, implies that all these 



11 



parameters can in principle be determined from the detailed analysis of the momentum 
dependence of the two-particle correlation functions and from the deviation of the one- 
particle invariant momentum distributions from the purely thermal behavior. 

For the special case when the measured QSCF-s are found to be momentum-independent 
and the IMD is well described by a Boltzmann momentum distribution, the a = case is 
concluded. In this case, we may measure the parameters At, R* and T* which are constants. 
The measurable parameters R* and T* are determined by four model parameters, Rg, Tq, 
t and b. Thus the model parameters cannot be determined uniquely from the IMD and 
QSCF observables in the a = special case. Inequalities can be obtained, along the lines of 
Ref. [Q], to restrict the values of the model parameters in this special case. For example, we 
have T < T* and Rg > R*- Note that a similar degeneracy of T* and R* has been found 
in Ref. [|IJ for the a = and 6 = 1 special case, where it was observed that the the same 
spectra and correlations can be obtained for different values of the three input parameters. 

The general results for the correlation function given in Eqs. ( |3B| - Po] ) indicate structural 
similarity with the Bose-Einstein correlation function parameters for a class of models which 
includes relativistic longitudinal flows, non-relativistic transverse flows and a transverse 
temperature profile fl,2~4}|. The radius parameters of the correlation function are momen- 
tum dependent both for the non-relativistic model presented here and for the model-class 
discussed in Refs. P,|2*3| . The structural similarity implies that the effective duration param- 
eter Atl , defined as the coefficient of in Eq. fl35|) , shall become momentum dependent: 
At* = At*(K) if the cooling effects are switched on. This new, momentum-dependent du- 
ration parameter At*(K) shall replace the momentum-independent At duration parameter 
in the correlation function if the temporal changes of the temperature are significant. 



V. APPLICATION TO NEUTRON AND PROTON INTERFEROMETRY 

In section [TV] it was discussed that the effective radius parameter R* (K) is not constant 
but depends on the momentum of the particles. Such an effect has indeed been seen in 
the measured proton-proton correlation functions in the 27 A1 ( 14 N, pp) reactions at E = 75 
MeV/nucleon 0,0: the larger the momentum of the protons the smaller the effective source 
size @. This feature is in qualitative agreement with the analytic results given in section 
PTI] , since the thermal radius Rt(K) is a decreasing function of the absolute value of the 
mean momentum, \K\. The effective source size is dominated by the smaller of the thermal 
and geometrical radius, thus it is also a decreasing function of the mean momentum of the 
particle pair. 

The results presented in the previous sections contain the essential ingredients which are 
needed to obtain a momentum-dependent radius parameter in the non-relativistic domain. 
However, to be utilized for nucleons emitted in intermediate energy heavy ion reactions 
|6|J7|], also two other important effects need to be considered, namely cooling of the source 
(see section [II B|) and final state interactions. Final state interactions can be included for 



pp pairs utilizing the Coulomb + strong interactions and Fermi-Dirac statistics, the strong 
interactions and Fermi-Dirac statistics for nn pairs and the strong interactions only for the np 
pairs, utilizing the Wigner-function formalism developed by S. Pratt and collaborators (see 
e.g. Refs. [0,|12j for the description of the method). Essentially, this calculation includes the 
evaluation of the two-particle relative wave-function with the above Coulomb and/or strong 
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interaction using the appropriate quantum statistical (anti)symmetrization and averaging 
the result over the particle distribution calculated from the model in section [II A . A Reid 
soft core potential is used to take into account the strong final state interactions. The 
essential consequences of the final-state interactions are that they completely modify the 
short-range parts of the two-particle correlations, and they create a peak at low relative 
momentum in the correlation functions. In the case of proton pairs, this peak is suppressed 
by the Coulomb repulsion which creates a hole centered at zero relative momentum in the 
pp correlation function. 

Light particle interferometry at intermediate energy heavy-ion collisions has been exten- 
sively studied, both experimentally and theoretically. Many attempts, using different models 
have been put forward, for example: simple Gaussian source parameterizations, evaporation, 
and transport models such as QMD and BUU (see Refs. and references therein). These 
models, containing different physical ingredients and information, have to different degrees 
been successful in describing the data. 

One of the purposes of this paper is to test if the model, presented in section |J, is 
applicable to nn and pp interferometry at intermediate energies. The model contains a rather 
small set of parameters, and by making simultaneous fits to n and p single energy spectra, 
and nn and pp correlation functions, rather hard constraints are put on the parameter set. 
Thus a qualitative and quantitative understanding of certain aspects of the source can be 
extracted. 

The evolution of the particle emission in a heavy-ion collision at intermediate energies 
may roughly be described as: production of pre-equilibrium particles; expansion and possible 
freeze-out of a compound source; possible evaporation from an excited residue of the source. 
Note though, that this separation is not very distinct and there is an overlap between the 
different stages. The importance of the various stages above also depends on the beam energy 
and the impact parameter of the collision. The model presented in this paper describes well 
the second stage above and, for long emission times, also part of the third stage. If these 
stages give the main contribution to the emission of nucleons, a satisfactory description of 
experimental data can be obtained. 





Rg (fin) 


T (MeV) 


a /t (fm/c)' 1 


b/to (fm/c)" 1 


d (fm/c)" 2 


Neutrons 


4.0 


3.0 


0.0 


0.018 


5.0 • 10~ 5 


Protons 


4.0 


5.0 


0.14 


0.036 


5.0 • 10~ 5 



TABLE I. Parameter values used for calculating n and p spectra and correlation functions 



We have applied our model to the reaction 40 Ar + 197 Au at 30 MeV/nucleon, to compare 
with experimental single spectra and correlation functions from the np correlation exper- 
iments described in Refs. J3 , 26| , 27 ] . With the parameter set presented in table | we have 
obtained a simultaneous fit to n and p single spectra, as well as nn and pp correlation 
functions. Thus we have used the same parameters for neutrons as for protons, except for 
the parameters To, a and b. As discussed in sections [I] and [Tl] the proton energy spectrum 
deviates from the thermal spectrum although the neutron spectrum is well approxi- 

mated by a Boltzmann distribution. The Coulomb interaction makes a difference between 
protons and neutrons, and this behavior is effectively obtained within our framework for 
protons by allowing a temperature gradient inside the source and a different flow parameter. 
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Protons emitted from the middle of the source roll down from a higher Coulomb-potential 
than those from the surface, thus protons from the middle are emitted with a higher kinetic 
energy than those emitted from the surface. This qualitative feature is similar to a system 
which is hotter in the middle than at the surface. Incorporating a temperature gradient 
inside the source does effectively describe such an effect. 

To keep the description as simple as possible in this first attempt to test the applica- 
bility of the model at intermediate energies, we have ignored impact parameter averaging 
and Coulomb interaction with the source. An averaging over different impact parameters 
could for example be done by allowing the radius parameter, Rq, to vary with the geomet- 
rical overlap of the projectile and target. However, such a prescription also introduces an 
additional uncertainty, and has therefore been ignored to keep the description simple. Thus 
the extracted source size will reflect a "mean" source size (assumed to be averaged over the 
different impact parameters). 

When performing the calculations we have taken into account the experimental energy 
thresholds and the acceptance region of the experimental set up. 



• Experiment (preliminary) 
□ Numerica 




Kinetic energy E k (MeV) 

FIG. 3. Experimental and numerically calculated proton energy spectrum. The filled circles 
represent preliminary inclusive proton emission data Q from the np correlation experiment de- 
scribed in Ref. [p6| 1, the open squares the full numerical simulation, while analytic approximations 
are represented by a solid curve (Eq. (|A2[)), dashed curve (Eq. (Xl)) and a dotted curve (Eq. (p6[)), 
respectively. The experimental data are preliminary, without fully estimated errors, why we have 
estimated the errors from the fluctuations of the data points. The parameter set in table | was 
employed for the numerical calculations. The energy distributions, f{Ek) are re-scaled so that 
f(E kfi ) = 1 for E kfi = 10 MeV. 



Experimental inclusive proton emission data (from the np correlation experiment 
described in Ref. p6fl ) are presented in Fig. [5] (filled circles) together with our fit (open 
squares). In addition we show the analytic approximate expression in Eqs. (p6|), ( |A1[) and 
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as a dotted, dashed and solid curve respectively. The neutron spectrum, not shown, 
is purely thermal with an effective temperature of 8.5 ± 1 MeV and is also well reproduced 
by our calculations. 

Proton correlation functions are presented in Fig. |as a function of the relative momen- 
tum q — \pi — pi\/2, integrated over the total pair momentum (p± +P2)- The experimental 
results are reproduced qualitatively, in some regions of the momentum space even quan- 
titatively. Considering an experimental uncertainty in q of 3.0 < Aq < 5.5 MeV/c P? 



which is not taken into account in our calculations, and the known difficulty to describe 
both the single-particle spectrum and the correlations for protons and neutrons at the same 
time ||10|| , we think that our result is interesting. The discrepancy for q < 10-15 MeV/c is 
partly due to the experimental uncertainty Aq. However, a discrepancy has also been seen 
in other models (see e.g. Ref. ||) and suggested explanations are the following: more than 
one source, Coulomb interaction with the source or pre-equilibrium emission. 



1 .2 - 



o 

c 
o 

o 
c 

c 

o 

_a 

CD 

o 
o 



0.8 - 



0.6 - 



0.4 - 



0.2 - 







0+0 ^ ni + 1 J 1 

a □'□ J a ~n at 



• Experiment 
□ Numerica 



10 20 30 40 50 60 

Relative momentum q (MeV/c) 

FIG. 4. Experimental and numerically calculated pp correlation function. The filled circles 
represent data from the np correlation experiment described in Ref. pq| , the open squares the full 
numerical simulation. The parameter set in table | was employed for the numerical calculations. 

The experimental nn correlation function contains some uncertainties, namely cross talk 
between the detectors, and experimental cuts at low relative momentum. From simulations, 
the cross talk (see Fig. 3e of Ref. ) has been found to mainly contribute to the "bump" 
around 15 < q < 40 MeV/c. Furthermore the sharp rise in C(q) at q ~ 10 MeV/c is 
based on rather few events with neutron kinetic energy smaller than about 10 MeV (see Fig. 
3d of Ref. [25|) close to detector energy threshold and acceptance limits. The momentum 
uncertainty, Aq, is smaller than 2 MeV/c. Considering these uncertainties the experimental 
C (q) is reproduced in an acceptable manner in Fig. |^. 

When performing the parameter fit we have found that it is not possible to simultane- 
ously reproduce n and p energy spectra, as well as nn and pp correlation functions, unless 
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10 20 30 40 50 60 

Relative momentum q (MeV/c) 

FIG. 5. Experimental and numerically calculated nn correlation function. The filled circles 
represent data from the np correlation experiment described in Ref. the open squares the full 
numerical simulation. The parameter set in table ffl was employed for the numerical calculations. 



the duration of the particle emission is taken large (several hundreds of fm/c). The long 
duration is found to be responsible for the low relative momentum behavior of the two- 
particle correlations. Thus we have used Hr(t) in Eq. flT%D for the time distribution. With 
the parameter set in table | we obtain numerically the mean emission time (t) and duration 
a(t): 

(t) « 520 fm/c cr(t) « 320 fm/c . 

(Note that (t) necessarily becomes large because of the broad time distribution limited to 
t > 0.) This scenario is in agreement with the general view of nuclear reactions in the 
energy domain of a few tenths of MeV/nucleon, namely a rather long-lived excited source 
emitting particles. Note though that a pure evaporative source gives a worse description 
of the experimental nn and pp correlation functions (see Ref. J^J) and fails to reproduce 
the single particle energy distribution. Thus the inclusion of the expansion seems to be 
mandatory, although the extracted flow velocities are small. 

The extracted geometrical source radius Rq ~ 4 fm is quite reasonable considering that a 
Gaussian parameterization is used (instead of a Woods-Saxon distribution, usually employed 
for the ground state density distribution). Note also that the extracted source size reflects 
an average over different times and impact parameters. The effective radius, R*, seen in the 
correlation function is smaller than Rq', 1.1 - 1.8 fm for protons (depending on momentum) 
and 2.5 fm for neutrons. 

From the experimental energy spectra we know that the effective temperature for neu- 
trons should be 8.5 ± 1 MeV and for protons in the range 15 - 25 MeV depending on the 
energy. The parameter set in table | gives 
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T cS = T* « 7.9 MeV 

for neutrons, and for protons, using the Taylor approximation in Eq. ( |A2|) , we obtain 

T cS « 20.1 MeV, £ fc = 10 MeV; T eff « 17.4 MeV, E k = 100 MeV . 

This is obtained with an input temperature of T = 3 MeV for neutrons, and for protons 
with a temperature profile that decreases from 5 MeV at the center of the source to 4.3 MeV 
at r = R G and 3.1 MeV at r = 2R G . 

The parameter set also gives the flow velocities 

u{n) « 0.018rc; u(p) « 0.036rc . 

The results presented in Figs. |] - 1] are obtained with cooling included according to Eq. 
17]. We have numerically found that the results change only to a minor extent if the cooling 



is excluded. Thus the different approximations in Eqs. (|2q), ( [All ) an d ( |A2j ) are useful also 
when cooling is included. 

A little note should be made on the parameter fit in our model. It is not possible to 
determine the parameters uniquely from the experimental data, though they are strongly 
constrained. Further constraints could be imposed by also investigating the dependence of 
the correlation function on the total pair momentum, though the current data set does not 
allow such an investigation. 



VI. SUMMARY 

In summary, we have calculated the invariant momentum distribution and the quan- 
tum statistical correlation function in a non-relativistic model for spherically symmetric 
non-relativistic expansion with possible temperature gradient, phenomeno logically includ- 
ing cooling and particle evaporation. We have thus presented a generalization of the previous 
work Jj]] by introducing new parameters. These parameters control the strength of the flow, 
the temperature gradient, the strength of cooling and the duration of particle emission. 

We have observed structural similarity of this non-relativistic model to the relativistic 
expansion described in Refs. Pl|33J, namely the effective radius parameters of the two-particle 
correlation functions and the effective slope parameters of the single particle spectra became 
momentum dependent due to the interplay of the local thermal scales and the geometrical 
scales. On this level, the model presented here is very similar to the relativistic one. 

The main effects of the temperature gradient are that it introduces i) a momentum- 
dependent effective temperature which is decreasing for increasing momentum, resulting in a 
suppression at high momentum as compared to the Boltzmann distribution; ii) a momentum- 
dependent effective source size which decrease with increasing total momentum. These 
qualitative features have been seen in non-relativistic heavy-ion collisions. 

The model presented here, including final state interactions, has been applied to mea- 
sured correlation functions and preliminary neutron and proton energy spectra in the re- 
action 40 Ar + 197 Au at 30 MeV/nucleon |[26|| . Agreement with the experimental data is 
obtained only if the duration time of the particle emission is large. The obtained param- 
eter set reflects a moderately large system (Gaussian radius parameter Rq = 4.0 fm) at a 
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moderate temperature (To(n) = 3 MeV and Tb(p) = 5 MeV) and small flow. Note, however, 
that the agreement between the model and the data was obtained only if the flow effect is 
included, i.e. within this phenomenological picture the inclusion of some flow is needed. 
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APPENDIX A: DETAILS OF ANALYTIC APPROXIMATIONS 

We have examined the validity of the saddle point approximation, presented in section 
[TTT] , both analytically and by comparing with a numerically generated particle distribution. 
The saddle point approximation normally yields a good approximation for large values of 
the expansion parameter (here a). However, for some parameter sets the saddle point 
approximation is not a good approximation for large a values (and large momenta p). This 
is because for large a, the emission function can develop more than one maximum, and the 
second order expansion around one of the maxima yields a poor approximation. The range 
of values of a and p for which the saddle point approximation is valid, depends on the other 
parameters used. 

For the case of two well separated maxima the saddle point approximation can be im- 
proved by summing up the contribution from the two maxima. Here we give the expression 
for the momentum distribution obtained in this approximation 



JVi(p)«c fl (27rRl(p)) im 




exp 




1 + a 2 E k Rl;/T tl 
1 + o?E k R%/T Q tl 



1 + a 2 E k R^ /T tl 
1 + a 2 E k R%/T tl 



(Al) 



An approximate expression can also be obtained for small values of the parameter a by 
making a Taylor expansion of the emission function around a = 0. Within this approxima- 
tion we obtain (identifying the integrated result as the first terms in an expansion of the 
exponential function): 
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^2+2 p2 p4 
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(A2) 



7*0 4t i? T0 Totg 

In the above expressions we have used the notation 

R T0 = R T (p = 0), i?* = i?*(p = 0) and T, = T„(p = 0) . 

Comparing these analytic approximations with numerical calculated momentum distri- 
butions we have found that for most parameter sets a good agreement can be found with 
either of the approximations. 
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